Susceptibilities with multi-quark interactions in PNJL model 



Abhijit Bhattacharyy£0 and Paramita DetQ 
Department of Physics, University of Calcutta, 
92, A. P. C. Road, Kolkata - 700009, INDIA 

Anirban LahirQ and Rajarshi RajH 

Center for Astroparticle Physics & Space Science, 
Bose Institute, Block-EN, Sector-V, 
Salt Lake, Kolkata-700091, INDIA 

& 

Department of Physics, Bose Institute, 
93/1, A. P. C Road, 
Kolkata - 700009, INDIA 

We have investigated the fluctuations and the higher order susceptibilities of quark number, isospin 
number, electric charge and strangeness at vanishing chemical potential for 2+1 flavor Polyakov loop 
extended Nambu-Jona-Lasinio model. The calculations are performed for the bound effective po- 
tential in the quark sector requiring up to eight quark interaction terms. These have been contrasted 
to the lattice results which currently have somewhat heavier quarks in the light flavor sector. The 
results show sufficient qualitative agreement. For comparison we also present the results obtained 
with the conventional effective potential containing upto six quark interaction terms. 
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I. INTRODUCTION 

Confinement and chiral symmetry breaking are the most fundamental properties of strong 
interaction physics at low temperature and density, where the physics is mainly governed by 
the non-perturbative QCD. In principle, the deconfinement phase transition and chiral phase 
transition are defined in two extreme limits of current quark mass. Deconfinement phase tran- 
sition and its order parameter is well defined for infinite current quark mass and chiral phase 
transition is exact for zero quark mass. But in real world with a finite value of quark masses, 
the nature of these two phase transitions is an open question. Another important feature of the 
QCD phase diagram is the existence of the critical end point (CEP), where first order phase 
transition, from hadronic phase to quark-gluon-plasma (QGP) phase, ends. However, the ex- 
act location of CEP is still unknown. Investigation of these properties of strongly interacting 
matter are necessary to understand the various astrophysical and cosmological scenario. 

The experimental explorations to understand the properties of strongly interacting matter 
has been studied at Relativistic Heavy- ion-collider (RHIC), Brookhaven [l| where the heavy 
nuclei collide with each other at relativistic energies to form hot and dense strongly interacting 
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matter. More data are expected from LHC and FAIR in future. But to analyze the data from 
the experiment, we need a thorough understanding of the theory of strong interaction physics. 

Due to our limited knowledge of the non-perturbative physics, QCD, which is the theory of 
strong interaction, can not be used to study the phase transition picture. In this regard Lattice 
QCD (LQCD) provides the most direct approach to study QCD at high temperature [2TH5I]. 
However LQCD has its own restrictions due to the discretization of space-time. Furthermore, 
at finite chemical potential, LQCD faces the well known sign problem. 

Another approach to study low energy limit of QCD and the QCD phase transition is to 
use effective theories of QCD. Polyakov loop extended Nambu-Jona-Lasinio model (PNJL) is 
one of the successful approach which combines the confinement and chiral symmetry breaking 
properties in a simple formalism. There have been series of work to study the thermodynamic 
pro per ties of strongly interacting matter using PNJL model for both 2 flavor and 2+1 flavor 
|16l - [24| . These studies suggest that this model reproduces the zero density lattice data quite 
successfully. However, the vacuum of the NJL part of the PNJL model seems to to be unbound 
in a 2+1 flavor scenario. A plausible solution of this problem has been proposed by Osipov et.al 
using ei ght- quark interaction term [25H281 ] . Also the 2 flavor PNJL model have been studied 
in Ref. [29l |30| with eight-quark interaction. In our previous work, we developed 2+1 flavor 
PNJL model with eight-quark interaction terms in the Lagrangian with three-momentum cutoff 
regularisation j3lj . 

The thermodynamic aspect of the phase transition from hadronic phase to QGP phase can be 
understood properly if we study the thermodynamic variables like quark number susceptibility 
(QNS), isospin number susceptibility (INS), specific heat (Cy) and speed of sound (v s ) etc. 
Susceptibilities are related to fluctuations via the fluctuation-dissipation theorem. A measure 
of the intrinsic statistical fluctuations in a system close to thermal equilibrium is provided by 
the corresponding susceptibilities. At zero chemical potential, charge fluctuations are sensitive 
indicators of the transition from hadronic matter to QGP. Also the existence of the CEP can 
be signalled by the divergent fluctuations. For the small net baryon number, which can be 
met at different experiments, the transition from hadronic to QGP phase is continuous and the 
fluctuations are not expected to lead any singular behavior. Recently, the computations on the 
lattice have been performed for many of these susceptibilities at zero chemical potentials [32T - 
35!]. It was shown that at vanishing chemical potential the susceptibilities rise rapidly around 
the continuous crossover transition region. 

The study of higher order moments of fluctuations are also necessary to locate the transition 
point more accurately. In 2 flavor QCD, it has been shown that the quark number and isospin 
number fluctuations increase with temperature and their fourth moments start to show pro- 
nounced peaks in the transition region from low to high temperature 0, [3^ | . In fact the higher 
order coefficients become increasingly sensitive in the vicinity of phase transition. Fluctuations 
are computed with respect to the quark chemical potential in 2 flavor PNJL model with three- 
momentum cutoff regularisation [22|, [37l - l39| . Also QNS at finite density has been estimated 
in some works within 2 flavor PNJL model [io|. Recently the idea of the Taylor expansion 
in terms of chemical potential for PNJL model has been computed within the constraint that 
the net strange quark density is zero, which is the case in ultra relativistic heavy ion collision 



41 1 . There have been recent calculations towards the fluctuations in 2+1 flavor Lattice QCD 
42h44| and also in the 2+1 flavor PNJL model taking upto six-quark determinant interaction 
terms [1^.|45|. Similar calculations have been carried out in Polyakov loop coupled quark-meson 
(PQM) model [4^ - |49j and its renormalization group improved version [Hoj . 

In this paper we have investigated the susceptibilities within 2+1 flavor PNJL model frame- 
work with two different kinds of NJL interaction. In one case we take the conventional form 
which takes into account upto six-quark determinant interaction, which we will call as model A. 
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In the other case an extra eight-quark interaction term is added to make the effective potential 
bound, which will be hereafter denoted by model B. We have also studied the specific heat and 
the speed of sound. Specific heat is related to the event- by- event temperature fluctuations [5l[ 
and mean transverse momentum fluctuations [52j in heavy-ion reactions. These fluctuations 
show diverging behavior near the critical end point (CEP). The speed of sound determines the 
flow properties in heavy-ion reactions [H, [54| . 

Our paper is organized as follows: In Sec. II, the basic formalism of the PNJL model and also 
the calculation of fluctuations have been discussed. Various thermodynamic quantities such as 
specific heat, speed of sound are also defined in this section. In the next section we describe 
our results and compared with the recent lattice data and also with other models like Polyakov 
extended quark- meson model (PQM). In the last section we conclude. 



II. FORMALISM 



A. Thermodynamic Potential 



The PNJL model was formulated to study the chiral properties and the confinement physics 
of the QCD phase transition at finite temperature and density. In this model quark dynamics 
is studied with a background gauge field having only the temporal c omp onent. For a detailed 
review of the PNJL model with 2 flavor and 2+1 flavor see Ref. [H, [H, H3, [H, [H, H^-SH HE H3l ■ 
Osipov et.al pointed out that the effective potential of the NJL part is seems to be unbound 
in the conventional form and they have introduced the eight-quark interaction terms in NJL 
model to stabilize the vacuum [25l428| . Also the 2 flavor PNJL model have been studied in Ref. 
[29. 30] with eight-quark interaction. In our previous paper we developed the 2+1 PNJL model 
with eight-quark interactions within three-momentum cutoff scheme [3l| for the first time. We 
have reproduced the thermodynamic properties of QCD, calculated on Lattice, at zero baryon 
density quite satisfactorily. The thermodynamic potential for the multi-fermion interaction in 
the mean field approximation (MFA) of the PNJL model can be written as [3l| , 
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gs and go are the four-quark and six-quark coupling constants respectively and g\ and gi are 
the eight-quark coupling constants. Here tjf — (ipfipf) denotes chiral condensate of the quark 

with flavor / and Ej = ^Jp 2 + Mj is the single quasi-particle energy. Here, constituent mass 

Mf of flavor / is given by the self-consistent gap equation; 



M f 



m f - 2g s uf 



9D „ , 2 

—<Jf+i<Jf+2 - 2gia f {a y 



- 4g 2 <J 3 f 



4 



where /, / + 1 and f + 2 take the labels of flavor u, d and s in cyclic order. So, when / = u then 
f + 1 = d and / + 2 = s and so on. In the above integrals, the vacuum integral has a cutoff A 
whereas the medium dependent integrals have been extended to infinity. Polyakov loop being 
the normalized trace of the Wilson line L, which is an SU(3) matrix, should lie in the range 
0s$ <f> <1. But earlier studies of PNJL model [H, 0, HI] show that the $ becomes greater 
than 1 above 2Tc- To solve this problem one has to take a proper Jacobian of transformation 
from the matrix valued field L to the complex valued field $. This will constrain the value of 
$ within 1. Thus one has to modify the Polyakov loop potential by introducing Vandermondc 
(VdM) term. The necessity of the modification was reflected through an excellent agreement of 
flavour mixing effects in PNJL model with lattice data as shown in Ref. [23| . There are a few 
extra terms in <l) in this prescription as compared to other recent works [l^, These 
extra terms are put on the basis of the global Z(3) symmetry of the Polyakov loop potential 
and therefore th eir p resence seem to be quite natural though may not be absolutely necessary 
as shown by [ijj l37l |. The modified potential W can be expressed as , 

"-W = "J*£n. KHmm (2) 

where U(<&, $, T) is the Landau- Ginzburg type potential given by [l7j . 

__!>£>,. _| (# . + # > ) + fc (4.,' (3) 

with, 

b 2 (T) =a + 0l (^) + a 2 (^) 2 + a 3 ( ^) 3 , (4) 

and 63, 64 being constants. J( < &, $) in eqn. ([2]) is known as VdM determinant (23|, is given by, 

J[$, !>] = (27/24tt 2 )(1 - 6$# + 4($ 3 + l> 3 ) - 3($$) 2 ) 

K is a phenomenological constant. Polyakov loop $ and its charge conjugate $ is defined as, 

$ = (Tr c L)/N c , $ = (Tr c Lt)/N c 

The parameter To is taken as 190 MeV, whereas the lattice determines its value to be 270 MeV 
for pure gauge theory. The reason to take a lower value of To is to get the crossover temperature 
(T c ) consistent with the lattice data. In this work we have taken the parameter set obtained 
in our previous paper [3~i| . For fixing the parameters m s ,A,gs,gD,gi,g2 we have used the 
following physical conditions; 

= 138 MeV f w = 93 MeV m K = 494 MeV f K = 117 MeV m n = 480 MeV m n , = 957 MeV 

and m u is kept fixed at 5.5 MeV. The parameters are given in table U for both PNJL model A 
and PNJL model B. 

For the Polyakov loop potential we choose the parameters which reproduces the lattice data 
of pure gauge thermodynamics [55j ]. It was shown in Ref. [13 that, pure gauge Lattice QCD 
data of scaled pressure, entropy and energy density are reproduced extremely well in Polyakov 
loop model using the ansatz ([3]) and (jj) with parameters summarized below, 



a = 6.75, ai = -1.95, a 2 = 2.625, a 3 = -7.44, b 3 = 0.75, 64 = 7.5 
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(MeV" 8 ) 


(MeV" 8 ) 




(MeV) 


ModelA 


5.5 


134.758 


631.357 


3.664 


74.636 


0.0 


0.0 


0.13 


181 


ModelB 


5.5 


183.468 


637.720 


2.914 


75.968 


2.193 


-5.890 


0.06 


169 



TABLE I: Parameters and Tc for model A and model B Lagrangians. 



B. Taylor expansion of pressure 



The pressure of the strongly interacting matter can be written as, 

P( T , AV MQ; Ms) = -tt{T, fi q , Hq, us), (5) 

where T is the temperature, fj, q is the quark chemical potential, fiQ is the charge chemical 
potential and [is is the strangeness chemical potential. From the usual thermodynamic relations 
we can show that the first derivative of pressure with respect to [i q gives the quark number 
density and the second derivative is the quark number susceptibility (QNS). 

Our first job is to minimize the thermodynamic potential numerically with respect to the 
fields u u , ad, <J S , and <l. The values of the fields can then be used to evaluate the pressure 
using the equation ([5]). Then we can expand the scaled pressure at a given temperature in a 
Taylor series for the chemical potentials [i q , fiQ, fxs as, 
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The flavor chemical potentials /i u , /id, Ms are related to [i q , fiQ, /is by, 
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Mu = Mg + gMQ, Md = f*q- gMQ 5 Ms = M<? - gMQ - Ms (8) 

It should be mentioned that one can also choose the independent chemical potentials as fi q , fii, 
fj,s- Then eqn.((8]) becomes, 

Mh = M<? + M/, Md = Mg - M/j Ms = M<?-Ms (9) 

where, /!/ is the isospin chemical potential. Eq.([6]) is a general expression of Taylor expansion 
of pressure for different chemical potentials. Since in this paper we are only concerned with the 
diagonal terms of the expansion, we can write the above equation in simpler way as, 

KT ' MX) -^C-(T)& (10) 
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Where X is q, Q or / and S. Here we will use the expansion around fix = 0, where the odd 
terms vanish due to CP symmetry. In this work we evaluate the expansion coefficients up to 
eighth order. To obtain the Taylor coefficients, first the pressure is obtained as a function of 
fix for each value of T, then fitted to a polynomial about fix = 0. All orders of derivatives are 
then obtained from the coefficients of the polynomial extracted from the fit. For the stability 
of the fit we have checked the values of least squares. 



C. Specific heat and speed of sound 



We have studied the specific heat Cy, which is important to find the location of CEP, the 
speed of sound, which determines the flow properties in heavy-ion reactions. The energy density 
e is obtained from the thermodynamic potential as, 
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The specific heat is defined as the rate of change of energy density with temperature at constant 
volume, which is given by, 



Cv ~df 



v dT 2 



For a continuous phase transition near CEP, it is expected that Cy shows a diverging behavior, 
which will translate into highly enhanced transverse momentum fluctuations or highly sup- 
pressed temperature fluctuations. The square of speed of sound at constant entropy S is given 
by, 
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(12) 



Divergence in specific heat near CEP means vanishing speed of sound. 



III. RESULT 



We now present the coefficients of the Taylor expansion of pressure for 2+1 flavor PNJL 
model with model A and model B Lagrangians and make a comparative study between the 
quark number susceptibility (QNS), isospin number susceptibility (INS), charge and strangeness 
susceptibility and their higher order derivatives. We then compare the specific heat and the 
speed of sound for both Lagrangians. We also compare our results with the recent lattice data 
available for 2+1 flavor with N T = 6 [43| and also with that of Polyakov extended quark-meson 
model (PQM) [Hfil. 



A. Coefficients of the Taylor expansion 



The pressure is fitted to a polynomial in fix using the "gnuplot" [5?} program at different 
values of temperature. Here we consider to take maximum eighth order term in the polynomial 
in fix- We restrict our expansion range to fi q ~ 300 MeV above which the diquark physics is 
expected to become important. Also the pion condensation and kaon condensation takes place 
in NJL model for fij > 70 MeV and fis > 240 MeV respectively. So we restrict our range within 
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FIG. 1: (color online). Variation of C2, C4, C6 and eg with T/Tc, for fix = fiq for models A and B. The 
arrows on the right show the corresponding SB limits. The lattice are data taken from Ref. [43| . 



jUj < 70 MeV and fis < 200 MeV below Tc. However, above Tc, approximate restoration of 
chiral symmetry implies that the chiral condensates become almost zero. So above Tc, we have 
extended the range of fi\ and fis for the better fit of the coefficients. Near Tc the x 2 (which is 
same as the least square here) of the fit varies rapidly with the variation of range of fix over 
which the fit was done. So near Tc we have fitted the pressure for 1 MeV gap of temperature 
and the data points are spaced by 0.1 MeV of chemical potential for all temperature values. 
The least-squares of all the fits came out to be 10 -10 or less. 

Now we study the behavior of the coefficients C2, C4, c§ and cs for three sets of chemical 
potentials for model A and model B. In figure ([T]) we show the variation of C2, C4, eg and c& 
with T/Tc for fix = fJ> q for both models and lattice data. It can be seen that QNS (c!]!) shows an 
order parameter like behavior. At low temperature there are small differences between model 
A and model B and model A is much closer to the lattice data. At high temperature c\ for 
model A reaches almost 98% of its ideal gas value whereas for model B it reaches almost 99% 
of its ideal gas value. However lattice data for N T — 6 at high temperature reaches almost 
the Stefan-Boltzmann (SB) limit. The fourth order derivative c| can be thought of as the 
susceptibility of c\. The figure shows a peak near Tc- Near Tc the model B shows much higher 
peak than the model A and the peak of the eight-quark interaction is closer to the lattice data. 
At higher temperature both cases match very well with the lattice data. But near and above 
2Tc lattice value converges to the SB limit, however both cases of PNJL model are slightly 
away from the SB limit. Note that both c\ and c| have only fermionic contribution in the 
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FIG. 2: (color online). Variation of C2, C4, C6 and cs with T/Tc, for \ix = \iq for models A and B. The 
arrows on the right show the corresponding SB limits. The lattice data are taken from Ref. [431 ]. 



SB limit. Since, the coupling strength is large enough for T < 2.5Tc, a sufficient amount of 
interaction is present in the system. So, it is expected that c| will not converge exactly to the 
SB limit within T < 2.5Tc- The higher order coefficients Cg and c| show interesting behavior 
near Tc- Although at very low and high temperatures both of them converge to zero. Near Tc, 
Cg shows sharp peaks for both cases. However for model B the peak is much sharper. Similar 
behavior can be observed for c\, which shows more peaks near Tc- The reason behind the 
peaks near the transition temperature may be due to the increase in fluctuation near Tc- The 
sharper peaks of model B is probably due to the introduction of enhanced repulsive interaction 
through eight-quark term. This was also reflected through the increase of scaled pressure in 
case of model B Lagrangian over the model A [3l|. The number of peaks increases near Tc 
for higher order coefficients. In figure @ the variation of susceptibility and the higher order 
coefficients for the charge chemical potential is shown. The nature of all the coefficients are 
same as the quark chemical potential. At high temperature the fluctuation for model B is 
closer to the SB limit compared to the model A. However lattice data is slightly above the SB 
limit for . For the case of , our data (for both cases) show a better convergence towards 
SB limit, unlike c\. At low temperature the behavior of model A is closer to the lattice data 
compared to the model B. The quartic fluctuations show a peak near Tc- The peak for model 
B is sharper than model A and the plot for model B matches well with the lattice result. The 
higher order coefficients show similar behavior as the quark chemical potential case. 

The figure (j3j) shows the variation of susceptibility and the higher order coefficients for the 
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FIG. 3: (color online). Variation of C2, C4, c$ and cs with T/Tc, for /j,x = /is for models A and B. The 
arrows on the right show the corresponding SB limits. The lattice are data taken from Ref. [431 ]. 



strangeness chemical potential with T/Tc- The cf for both the models are slightly different 
from the lattice data. Both the plots are almost 98% of the SB limit at high temperature, 
however the lattice data coincides with the SB limit. The cf has a similar behavior as the c\. 
However the peak is not near Tc in this case. Near Tq we can see a small bump for both type 
of Lagrangian, but the peak in both case is at higher temperature. This is due to the fact that 
during the chiral crossover the strange quark (the only element which carries strangeness) is 
sufficiently heavy and the corresponding condensate a 8 melts at much higher temperature than 
Tq- The maxima of da s j dT and cf coincide at the temperature where peaks are observed. This 
behavior is quite consistent with the lattice results which also indicates two peaks. However 
in case of lattice, the peak at Tc is higher than the second peak. So, one can not really pin 
down the cause of the double peak structure. It may be a model artifact. At high temperature 
both the curves are above the SB limit. But the lattice data is slightly below the SB limit. 
The higher order derivatives also show a sharp peak at the transition temperature followed by 
a broader peak at higher temperature for both model A and Model B. 

For the sake of completeness we have also plotted different moments of pressure for isospin 
chemical potential /ij in fig. Q. It is clearly seen that, c\ has also an order parameter like 
behavior like all other susceptibilities both models. At high temperature, the plot for model B 
reaches almost 99.5% of the SB limit, whereas the plot for model A reaches close to 98% . For 
C4, the plot for the model B shows higher peak than the case of model A at Tc and both of 
them converges very well towards the SB limit at high temperature. The c\ and c\ shows very 
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FIG. 4: (color online). Variation of C2, C4, C6 and cs with T/Tc, for = fj,i for models A and B. The 
arrows on the right show the corresponding SB limits. The lattice data are taken from Ref. [431 ]. 



rapidly fluctuating peak structure near Tc and goes to zero in both high and low temperature 
regime. 

In figure ([5]) we show a comparative study for C2, C4, eg and eg, with T /Tc for (j, q , us and 
/xq for both type of Lagrangian. In all cases we can see that the quadratic fluctuations rise 
rapidly in the transition region where the quartic fluctuations show a peak. This peak is most 
pronounced in case of \i q . The generic form of this temperature dependence, a smooth crossover 
for quadratic fluctuations and a peak in quartic fluctuations, is in fact expected to occur in 
the vicinity of the chiral phase transition of QCD. For all of the coefficients the fluctuation 
in [i q direction is strongest followed by the isospin direction fluctuations and the strangeness 
fluctuation. The charge fluctuation is the least pronounced. The quartic fluctuation for fis 
shows a peak at much higher temperature than the transition temperature. However other 
three c|, C4 and c\ show peaks at Tq. For higher order coefficients the strangeness and charge 
fluctuations are negligible near Tc compared to the quark number and isospin fluctuations. 

In fig ^ we have plotted the kurtosis i.e. the ratio of c£ /c£ for both type of potential 
(where X = q,Q or S) and compared with the lattice data. The plot for c\jc\ for model B 
shows more fluctuation near Tc than model A. However lattice data shows higher fluctuation 
near Tc than the model study. At higher temperature both models coincide with the lattice 
data and converges well with the SB limit. Our results are qualitatively similar with PQM 
result [48]. In case of the ratio c^jc^ the model B shows more fluctuation than the model A 
and as well as the lattice data. The model B shows almost 99% convergence with the SB limit 
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FIG. 5: (color online). Comparative study of e 2 , cf , c§ and c$ for models A and B, where X=q 
I or S. 
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arrows on the right show the corresponding SB limits. The lattice data are taken from Ref. The 
upper left panel corresponds to the quark chemical potential, the upper right panel corresponds to the 
charge chemical potential and the lower panel corresponds to the strangeness chemical potential 



at high temperature. For the strangeness fluctuation, we can see two peaks in cf / cf curve for 
the both models. First peak occurs at chiral transition for light flavors and second peak occurs 
when chiral transition occurs in strange sector. At intermediate temperatures PNJL model 
overestimates the ratio than LQCD result. This feature is also observed in PQM model [Hj]. In 
PQM model the ratio approaches the SB limit at high temperature and near the transition region 
sharp peaks appear. In our case the ratio shows sharper peak near the transition temperature 
and a broader peak at a higher temperature than Tc- However lattice data shows much higher 
fluctuations than the plots for both models near Tc- At high temperature the model A is closer 
to the SB limit than the model B. 

B. Specific heat and the speed of sound 

We now discuss the thermodynamic quantities like specific heat (Cy) and the speed of sound 
(v s ). In Fig. ([7]) we have plotted CV/T 3 with T/Tc for both models. From the plot we can see 
that Cy grows with increasing temperature and reaches a peak at Tq for both models. However 
model B shows a sharper peak at Tc compared to the model A. Just above Tc both the plots 
decrease sharply for a short range of temperature. Thereafter it gradually converges to a value 
slightly lower than the ideal gas value at high temperature. However the convergence towards 
the ideal gas value is better in case of model B. For comparison we have also plotted the values 
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FIG. 7: (color online). Variation of Cv/T 3 , ie/T 4 , v 2 s and p/e with T/T c , for models A and B. The 
arrows on the right show the corresponding SB limits. The upper left panel shows the plot for Cv/T 3 
and 4e/T 4 with T/Tg, the upper right panel shows the comparison of vl for models A and B and the 
lattice data given by [4411 and the bottom panel shows the comparison of p/e of both models and the 
lattice data given by [42l . |44| 



of 4e/T 4 , at which the specific heat is expected to coincide for a conformal gas. From the graph 
we can see that the specific heat converges very well with the 4e/T 4 at high temperature for 
both the models. Both the plots of Cv/T 3 and 4e/T 4 show similar behavior as PQM model 
[48j | . In PQM model Cy grows with temperature and shows a sharp peak at Tq- After the 
peak a broad bump is found around 1.2Tc an d then Cy goes gradually to the ideal gas value. 

We now consider the speed of sound and p/e for model A and model B in figure (J7J. We can 
see that v 2 s is slightly below the ideal gas value at temperature 2.5Tc for both cases. Our result 
is quite consistent with the lattice data for 2+1 flavor staggered fermions reported in (5(| . A 
similar behavior as ours is observed in PQM model (48[. We get the minimum of «f just below 
the Tc similar to lattice data [56| and the softest point of the equation of state is found to be 
{p/e) m in ~ 0.07 for model A and {p/t) m in ~ 0.06 for model B. Model B gives better agreement 
with lattice data, which has its softest point of equation of state as (p/e) m in ~ 0.05. We have 
compared two sets of lattice data of Ref. l42l a nd Ref. [H| with our model study. Our result 
shows a better agreement with that of Ref. |42| . The softest point of equation of state of [42[ 
is at (p/e) m i n w 0.08, whereas the softest point of Ref. [44[ is at much higher value ~ 0.13. In 
case of PQM model the softest point of the equation of state is found to be around « 0.04 [48[ 
which is closer to our value. 
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IV. DISCUSSION 



We have studied the various fluctuations and some of the thermodynamic quantities using 
two diferent versions of PNJL model to understand the properties of the strongly interacting 
matter. In fact it is expected that the susceptibilities and the higher order fluctuations might 
provide the direct evidence of the order of the QCD phase transition. A pronounced peak in 
the susceptibility can depict the crossover transition and the sharp diverging behavior would 
indicate the existence of a phase transition. 

We have also obtained the susceptibilities and the higher order derivatives by the Taylor 
expansion of pressure for two kinds of PNJL model near /ix — 0, where X = q, I or Q and 
S. In all cases the second derivative of pressure which is known as the susceptibility, show a 
steep rise near the transition region, which indicates near the transition region the fluctuation 
increases. However at higher temperature c£ almost saturates and almost converges to the ideal 
gas value. This result is quite consistent with the lattice data. The higher order fluctuation c 4 
shows a peak near Tq for both models and the result matches with the lattice data. The finite 
height of the peak confirms the crossover nature of transition at /x = 0. Both c 6 and c 8 show 
rapid variation around Tq for all cases. 

We have also calculated the specific heat, speed of sound for both kinds of PNJL model. The 
plot for Cy/T 3 , after showing a peak at Tc, converges very well to 4e/T 4 curve at high T. At 
high temperature v"l almost reaches its ideal gas value 1/3 and the softest point of the equation 
of state has a better agreement with lattice result for model B where eight-quark interaction is 
taken into account. 

In our formalism, we do not include pion condensate, kaon condensate and diquark condensate 
which may play an important role for higher values of chemical potential. Inclusion of those 
degrees of freedom may improve our result. Involved though, such studies are in progress. 
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